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Specification 
Title of the Invention 
Limit Cycle Autotuning Method and 
Heat /Cool Control Apparatus 

5 

Background of the Invention 

The present invention relates to a heat /cool 
control technique of performing temperature control by 
properly switching a heat mode of outputting a 

10 manipulated variable to a heating actuator and a cool 
mode of outputting a manipulated variable to a cooling 
actuator and, more particularly, to a limit cycle 
autotuning method and heat /cool control apparatus which 
adjust control parameters by generating a limit cycle 

15 with a constant manipulated variable amplitude. 

A control technique called heat/cool control 
which selectively uses the heating capacity of a heater 
and the cooling capacity of a cooler is used for 
temperature control on air in the chamber of a 

20 thermostatic chamber unit, temperature control on a 

plastic material in an injection molding machine, and 
the like. The typical arrangement of an apparatus using 
this technique includes a mechanical unit which can heat 
and cool air and a fan which sends air with an arbitrary 

25 temperature at a proper airflow and sends heated or 
cooled air into a chamber/room as a control target . 

Fig. 9 shows the arrangement of such a 
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heat /cool control system. Reference numeral 21 denotes 
a thermostatic chamber; 22, a heater which heats air; 
23, a cooler which cools air; 24, a temperature sensor 
which measures the temperature in the chamber; and 25, a 
5 heat/cool control apparatus. The heat/cool control 

apparatus 25 computes a manipulated variable MV on the 
basis of a set point ( intra- chamber temperature set 
point) SP set by an operator and a controlled variable 
(intra -chamber -temperature) PV measured by the 

10 temperature sensor 24, and outputs the result to the 
heater 22 or cooler 23. The heating capacity of the 
heater 22 is increased/decreased by an SCR in accordance 
with the manipulated variable MV. The cooling capacity 
of the cooler 23 is increased/decreased by an inverter 

15 in accordance with the manipulated variable MV. 

In brief, a control technique for the 
heat /cool control system in Fig. 9 is a temperature 
control technique which implements one PID control 
system and is of a heating/cooling switching type that 

20 operates the heater 22 in accordance with the 

manipulated variable MV (heat mode) if the manipulated 
variable MV is 0% or higher and operates the cooler 23 
in accordance with the manipulated variable MV (cool 
mode) if the manipulated variable MV is lower than 0%, 

25 as shown in Fig. 10. As described above, the heat /cool 
control is a control technique which implements a single 
loop control system which is applied to an apparatus 
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arrangement including heating and cooling actuators in 
one loop temperature control system and handles the 
heating and cooling actuators as a set. According to 
this technique, the manipulated variable MV is output to 
5 a heat -side operation end (heating actuator) in the heat 
mode, whereas the manipulated variable MV is output to a 
cool- side operation end (cooling actuator) in the cool 
mode. 

In the above heat /cool control system, the PID 

10 parameters in the heat/cool control apparatus 25 must be 
switched to proper values in each of the heat mode and 
the cool mode. In order to automatically obtain this 
proper value, an autotuning method has been proposed. 
This autotuning method is disclosed in, for example, 

15 Japanese Patent Laid-Open No. 5-289704. The autotuning 
method uses a limit cycle scheme. First of all, PID 
parameters on the heat mode side are obtained by the 
general limit cycle scheme of operating only the heater 
(in the first half of autotuning). The heater and 

20 cooler are then operated to generate a limit cycle (in 
the second half of autotuning) . The proportional band 
of PID parameters on the cool mode side is calculated on 
the basis of a limit cycle amplitude in the first half 
of autotuning and a limit cycle amplitude in the second 

25 half of autotuning. Thereafter, autotuning is 
terminated. 

Some apparatus using heat /cool control is 
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designed such that forced heating by a heater is 
required to raise the temperature, and forced cooling by 
a cooler is required to lower the temperature. Some 
thermostatic test apparatus is used at room temperature 
5 (about 25°C) . In such an apparatus, autotuning is 

performed near room temperature. Assume that a limit 
cycle is to be generated by using only the heater at 
this time. In this case, after the intra-chamber 
temperature is raised, the temperature must be lowered 

10 while the apparatus is left standing at room 

temperature. As a consequence, a condition occurs in 
which the temperature hardly drops, i.e., a condition 
occurs in which it is impossible to generate a limit 
cycle. Likewise, assume that a limit cycle is to be 

15 generated by using only the cooler. In this case, after 
the intra-chamber temperature is lowered, the 
temperature must be raised while the apparatus is left 
standing at room temperature. As a consequence, a 
condition occurs in which the temperature hardly rises. 

20 If a limit cycle cannot be generated by using 

only the heater or cooler, since it is difficult to 
generate a limit cycle by operating only the heater as 
in the above autotuning method, the conventional 
autotuning method cannot be used. Conventionally, as 

25 described above, in the heat/cool control system which 
cannot generate a limit cycle by using only the heater 
or cooler, it is impossible to obtain control parameters 
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suitable for the heat mode and cool mode, respectively. 

The above description is based on the premise 
that autotuning is performed near room temperature. 
This is, however, an example of autotuning. In some 
5 case, when a chamber is left standing after the 

intra- temperature is raised, the equilibrium point of 
temperature is settled at a temperature higher than room 
temperature due to a surrounding environment or the heat 
capacity of a furnace body. In such a case, in order to 

10 lower the temperature, forced cooling by a cooler is 
required. In contrast to this, in some case, when a 
chamber is left standing after the intra- temperature is 
lowered, the equilibrium point of temperature is settled 
at a low temperature. In such a case, in order to raise 

15 the temperature, forced hating by a heater is required. 
Therefore, when autotuning is to be performed at a 
temperature other than room temperature, the above 
problem can arise. 
Summary of the Invention 

20 The present invention has been made to solve 

the above problem, and has as its object to provide a 
limit cycle autotuning method and heat/cool control 
apparatus which can obtain control parameters suitable 
for both the heat mode and the cool mode even if an 

25 autotuning function of heat /cool control cannot generate 
any limit cycle by using a heater or cooler alone. 

In order to achieve the above object. 
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according to the present invention, there is provided a 
limit cycle autotuning method of calculating a control 
parameter by alternately performing operation of 
outputting a predetermined heat -side manipulated 
5 variable set point to a heating actuator and operation 
of outputting a ^predetermined cool -side manipulated 
variable set point to a cooling actuator in a heat/cool 
control apparatus which performs temperature control by 
performing feedback control computation with respect to 

10 a deviation between a set point and a controlled 

variable on the basis of the control parameter, and 
properly switching a heat mode of outputting a 
manipulated variable to the heating actuator and a cool 
mode of outputting a manipulated variable to a cooling 

15 actuator, comprising the first limit cycle generation 
step of generating a first limit cycle of alternately 
outputting the heat -side manipulated variable set point 
and the cool-side manipulated variable set point, the 
first control response detection step of detecting a 

20 first control response corresponding to the first limit 
cycle, the second limit cycle generation step of 
generating a second limit cycle by changing one of the 
heat -side manipulated variable set point and the 
cool- side manipulated variable set point on the basis of 

25 predetermined change instruction information for 

instructing which one of the heat -side manipulated 
variable set point and the cool- side manipulated 
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variable set point is to be changed after the first 
limit cycle and a predetermined manipulated variable 
change ratio indicating a degree of the change, the 
second control response detection step of detecting a 
5 second control response corresponding to the second 

limit cycle, and the control parameter calculation step 
of calculating the control parameter for each of the 
heat mode and the cool mode on the basis of the detected 
first and second control responses. 
10 Brief Description of the Drawings 

Fig. 1 is a block diagram showing the 
arrangement of a heat /cool control apparatus according 
to the first embodiment of the present invention; 

Fig. 2 is a flow chart showing the operation 
15 of the heat /cool control apparatus in Fig. 1 at the time 
of execution of limit cycle autotuning; 

Fig. 3 is a flow chart showing the operation 
of the heat /cool control apparatus in Fig. 1 at the time 
of execution of limit cycle autotuning; 
20 Fig. 4 is a flow chart showing the first 

ascent/descent extreme value detection processing by a 
limit cycle autotuning computation section in Fig. 1; 

Fig. 5 is a timing chart for explaining the 
first ascent/descent extreme value detection processing; 
25 Fig. 6 is a timing chart for explaining 

controlled variable amplitude calculation processing; 

Fig. 7 is a block diagram showing the 
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arrangement of a heat /cool control apparatus according 
to the second embodiment of the present invention; 

Fig. 8 is a flow chart showing the operation 
of the heat/cool control apparatus in Fig. 7 at the time 
5 of execution of a limit cycle autotuning; 

Fig. 9 is a block diagram showing the 
arrangement of a conventional heat/cool control system; 
and 

Fig. 10 is a graph for explaining manipulated 

10 variable outputs in heat/cool control. 

Description of the Preferred Embodiments 

The present invention will be described in 
detail next with reference to the accompanying drawings. 
First Embodiment 

15 The principle of limit cycle autotuning which 

automatically adjusts PID parameters according to a 
limit cycle with a constant manipulated variable 
amplitude will be described first. By continuing a 
manipulated variable output equivalent to ON-OFF 

20 control, a ascent /descent phenomenon (limit cycle) of a 
controlled variable PV corresponding to the 
characteristics of a control target can be obtained. In 
this case, a transfer function P of the control target 
is given by 

25 P = Kpexp(-Lps)/(1 + Tps) ...(1) 

where Kp is a process gain, Tp is a process time 
constant, Lp is a process dead time, and s is a Laplace 
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operator. 

Let AU be the upper limit value of a 
manipulated variable MV with a constant amplitude which 
is output at the time of execution of limit cycle 
5 autotuning, and AL be the lower limit value of the 

manipulated variable. If the controlled variable PV is 
equal to or lower than a set point SP, manipulated 
variable MV = AU is output to the control target. If 
the controlled variable PV is higher than the set point 

10 SP, manipulated variable MV = AL is output. If a limit 
cycle is generated when the controlled variable PV sets 
the time average of manipulated variable outputs in the 
limit cycle to almost 0.5 (AU + AL), a maximum change 
ratio dPVmax of the controlled variable PV becomes 

15 0 . 5 ( Kp/Tp ) ( AU - AL). 

A necessary time (manipulated variable 
switching elapsed time) Tha from the instant at which 
the manipulated variable MV is switched from AL to AU to 
the instant at which the controlled variable PV reaches 

20 the lowest descent point is equal to the process dead 
time Lp. Likewise, a necessary time (manipulated 
variable switching elapsed time) Thb from the instant at 
which the manipulated variable MV is switched from AU to 
AL and the instant at which the controlled variable PV 

25 reaches the highest ascent point is also equal to the 
process dead time Lp. When the time average of 
manipulated variable outputs in the limit cycle becomes 
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0.5(AU + AL), the necessary time from the lowest descent 
point of the controlled variable PV to the highest 
ascent point becomes almost 2Lp. According to the above 
description, an amplitude Ah of the controlled variable 
5 PV can be approximately given by 

Ah = 2LpdPVmax = (Kp/Tp)Lp(AU - AL) ...(2) 
Since an ordinary control target has a slight 
higher-order lag characteristic, the actual amplitude Ah 
is about 75% of the value given by equation (2) as 
10 indicated by the following equation: 

Ah = 0.75(Kp/Tp)Lp(Au - AL) ...(3) 
In this case, as guidelines for setting PID 
parameters including a proportional band Pb, integral 
time Ti, and derivative time Td, the following are 
15 generally proposed: 

Pb - <5KpLp/(aTp) .-.(4) 
Ti = j3 Lp . . . ( 5 ) 

Td = T Lp . . . ( 6 ) 

In equations (4) to (6), the constants a, £ , r, and 6 
20 are, for example, 1.2, 2, 0.42, and 100, respectively. 

When a limit cycle is generated, the amplitude 
Ah of the controlled variable PV and the manipulated 
variable switching elapsed times Tha and Thb (= Lp) are 
obtained. Letting Th be the average of Tha and Thb or 
25 an appropriate one of them and substituting them into 
equations (4) to (6), PID parameters are calculated as 
follows : 
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Pb = <5KpLp/(aTp) = <5Ah/{£(AU - AL) } 
Ti = )3Lp = j3Th 



• . . (8) 



. • .(7) 



Td = 7Lp = TTh 



...(9) 



In equation (7), the constant S is 0.9. 
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The principle of a heat/cool control apparatus 



according to this embodiment will be described next. 
The value of the manipulated variable MV to be output to 
the heater during a limit cycle is set to a heat -side 
manipulated variable set point ATh (0 < ATh ^ 100%), 

10 and the value of the manipulated variable MV to be 

output to the cooler is set to a cool-side manipulated 
variable set point ATc (-100 < ATc ^ 0%). First of 
all, in the first limit cycle, if the controlled 
variable PV is equal .to or lower than the set point SP, 

15 manipulated variable MV = ATh is output to the heater. 
If the controlled variable PV is higher than the set 
point SP, manipulated variable MV = ATc is output to the 
cooler. In the second limit cycle, the heat -side 
manipulated variable set point is changed to MATh (0 < 

20 U < 1). 



waveform (controlled variable PV) in the first limit 
cycle differs from that in the second limit cycle. This 
amplitude difference is related to the process gain 
25 (heat -side process gain) set when the heat -side 

manipulated variable set points ATh and M ATh are output 
and the process gain (cool-side process gain) set when 



In this case, the amplitude of a response 
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the cool- side manipulated variable set point ATc is 
output . 

The ratio between a heat -side process gain and 
a cool- side process gain is also the ratio between a 
5 proportional band (heat -side proportional band) Pb_Heat 
suitable for the calculation of the manipulated variable 
MV to be output to the heater and a proportional band 
(cool-side proportional band) Pb_Cool suitable for the 
calculation of the manipulated variable MV to be output 
10 to the cooler. This ratio K is determined as follows: 
Pb_Heat : Pb_Cool = K : 1 - K ...(10) 
According to the ratio K , an amplitude Al of 
the controlled variable PV in the first limit cycle and 
an amplitude A2 of the controlled variable PV in the 
15 second limit cycle are given by 

Al = A { K ATh + (1 - /C)ATc} ...(11) 
A2 = A{/Z/CATh + (1 - K )ATc} ...(12) 
In equations (11) and (12), A is a proportional 
coefficient . 

20 If the proportional coefficients A are 

eliminated from equations (11) and (12), and the 
equations are combined with respect to the ratio k , the 
following is obtained: 

K = {ATC(A2 - Al)} 
25 /{ATc(A2 - Al) + ATh(A2 - MAI)} ...(13) 

If a cool- side manipulated variable set point 
M ATc (0 < ii < 1) is changed instead of the heat-side 
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manipulated variable set point ATh in the second limit 
cycle, the following is obtained: 
K = {ATc( JJLK1 - A2)> 

/{ATC(A1 - A2) + ATc(MAl - A2)> ...(14) 
5 In this case, if an average proportional band 

Pb_Ave which is the average of the heat -side 
proportional band Pb_Heat and cool -side proportional 
band Pb_Cool is simply calculated on the basis of the 
amplitude A2 of the controlled variable PV which appears 
10 when the heat-side manipulated variable set point is 
changed to M ATh in the second limit cycle, the 
following is obtained: 

Pb_Ave = <5 A2/{ 6 ( MATh - ATc ) } ...(15) 
Likewise, if the average proportional band 
15 Pb__Ave is simply calculated on the basis of the 

amplitude A2 of the controlled variable PV which appears 
when the cool-side manipulated variable set point is 
changed to tt ATh in the second limit cycle, the 
following is obtained: 
20 Pb_Ave = <5A2/{6(ATh -/XATc)} ...(16) 

The relevant ratio of the heat -side process 
gain (heat -side proportional band) to the average 
proportional band Pb_Av is K , and the relevant ratio of 
the cool- side process gain (cool- side proportional band) 
25 is 1 - k . When, therefore, the heat -side manipulated 

variable set point which has been ATh in the first limit 
cycle is changed to tt ATh in the second limit cycle, or 
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the cool-side manipulated variable set point which has 
been ATc in the first limit cycle is changed to MATc in 
the second limit cycle, the heat -side proportional band 
Pb_Heat and cool- side proportional band Pb_Cool are 
5 obtained as follows : 

Pb_Heat = 2 K Pb_Ave ...(17) 
Pb_Cool = 2(1 - K)Pb_Ave ...(18) 
Let Th__Heat be a manipulated variable 
switching elapsed time from the instant at which the 

10 manipulated variable MV is switched from the cool- side 
manipulated variable set point ATc (or /XATc) to the 
heat -side manipulated variable set point MATh (or ATh) 
in the second limit cycle to the instant at which the 
controlled variable PV reaches the minimum value, and 

15 Th_Cool be a manipulated variable switching elapsed time 
from the instant at which the manipulated variable MV is 
switched from the heat -side manipulated variable set 
point MATh (or ATh) to the cool-side manipulated 
variable set point ATc ( tt ATc ) to the instant at which 

20 the controlled variable PV reaches the maximum value. 

In general heat/cool control process apparatus 
design, there are no great differences in integral and 
derivative times as PID parameters between the heat mode 
and the cool mode. When, therefore, either the 

25 heat -side manipulated variable set point in the second 
limit cycle is to be changed to l± ATh or the cool- side 
manipulated variable set point is changed to fi ATc , an 
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average integral time Ti_Ave and average derivative time 
Td_Ave are obtained from the average of the manipulated 
variable switching elapsed times Th_Heat and Th_Cool. 
The obtained average integral time Ti_Ave and average 
5 derivative time Td_Ave may be used as PID parameters on 
the heat mode side and cool mode side, respectively. 

That is, the average integral time Ti_Ave and 
average derivative time Td__Ave can be calculated as 
follows by applying the manipulated variable switching 
10 elapsed times Th_Heat and Th_Cool to equations (8) and 
(9): 

Ti_Ave = )3 (Th_Heat + Th_Cool)/2 ... (19) 

Ti_Ave = 7(Th_Heat + Th_Cool)/2 ...(20) 
As an integral time (heat -side integral time) 
15 Ti_Heat suitable for the calculation of the manipulated 
variable MV to be output to the heater and an integral 
time (cool- side integral time) Ti_Cool suitable for the 
calculation of the manipulated variable MV to be output 
to the cooler, the average integral time Ti_Ave may be 
20 simply used. As a derivative time (heat -side derivative 
time) Td_Heat suitable for the calculation of the 
manipulated variable MV to be output to the heater and a 
derivative time (cool-side derivative time) Td_Cool 
suitable for the calculation of the manipulated variable 
25 MV to be output to the cooler, the average derivative 
time Td_Ave may be simply used. 

Ti_Heat = Ti_Ave ...(21) 
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Ti_Cool = Ti_Ave ... (22) 

Td_Heat = Td_Ave ... (23) 

Td_Cool = Td_Ave ... (24) 

The arrangement of the heat/cool control 
5 apparatus according to this embodiment will be described 
on the basis of the above principle. Fig. 1 shows the 
arrangement of the heat /cool control apparatus according 
to the first embodiment of the present invention. The 
heat /cool control apparatus in Fig. 1 includes a set 

10 point input section 1 which inputs the set point SP set 
by the operator of the control apparatus , a controlled 
variable input section 2 which inputs the controlled 
variable PV detected by a sensor (not shown) , a 
manipulated variable output section 3 which outputs the 

15 manipulated variable MV to a heater (not shown) serving 
as a heating actuator for implementing a heating 
function for temperature control or a cooler (not shown) 
serving as a cooling actuator for implementing a cooling 
function for temperature control, a PID control 

20 computation section 4 which calculates the manipulated 
variable MV by performing PID control computation for 
the deviation between the set point SP and the 
controlled variable PV on the basis of control 
parameters (PID parameters), a heat-side manipulated 

25 variable storage section 5 which stores a predetermined 
heat -side manipulated variable set point ATh (0 < ATh ^ 
100%) to be output to the heater at the time of 
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execution of autotuning, and a cool- side manipulated 
variable storage section 6 which stores a predetermined 
cool-side manipulated variable set point ATc (-100 < ATc 
^ 0%) to be output to the cooler at the time of 
5 execution of autotuning. 

The heat/cool control apparatus further 
includes a limit cycle autotuning computation section 7 
which generates the first limit cycle in which the 
heat -side manipulated variable set point ATh and 

10 cool- side manipulated variable set point ATc are 

alternately output and the second limit cycle in which 
either the heat -side manipulated variable set point ATh 
or the cool- side manipulated variable set point ATc is 
changed, and detects the first control response 

15 corresponding to the first limit cycle and the second 
control response corresponding to the second limit 
cycle, a manipulated variable change ratio storage 
section 8 which stores predetermined change instruction 
information CI for instructing which one of the 

20 heat-side manipulated variable set point ATh and the 
cool -side manipulated variable set point ATc is to be 
changed after the first limit cycle, and a predetermined 
manipulated variable change ratio U (0 < fi < 1) 
indicating the degree of this change, a two-step limit 

25 cycle computation section 9 which instructs the limit 
cycle autotuning computation section 7 to execute the 
first and second limit cycles, a first limit cycle 
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information storage section 10 which stores the 
information of the first control response corresponding 
to the first limit cycle, a second limit cycle 
information storage section 11 which stores the 
5 information of the second control response corresponding 
to the second limit cycle, and a PID parameter 
calculation section 12 which calculates PID parameters 
for the heat and cool modes on the basis of the detected 
first and second control responses, and sets the 
10 calculated PID parameters in the PID control computation 
section 4. 

The PID control computation section 4 
implements a control computation means. The manipulated 
variable change ratio storage section 8 implements a 

15 manipulated variable change ratio storage means. The 
limit cycle autotuning computation section 7 and 
two-step limit cycle computation section 9 implement a 
limit cycle generating means and control response 
detection means, respectively. The PID parameter 

20 calculation section 12 implements a control parameter 
calculation means . 

The operation of the heat/cool control 
apparatus in Fig. 1 will be described with reference to 
Figs. 2 and 3. Referring to Figs. 2 and 3, steps 102 to 

25 105 show a sequence for generating the first limit 

cycle; steps 106 and 107, a sequence for detecting the 
first control response; steps 108 to 112 and 115 to 117, 



- 18 - 



a sequence for generating the second limit cycle; steps 
113, 114, 118, and 119, a sequence for detecting the 
second control response; and steps 120 to 122, a 
sequence for calculating control parameters. 
5 The set point SP is set by the operator of the 

autotuning apparatus and input to the PID control 
computation section 4, limit cycle autotuning 
computation section 7, and two-step limit cycle 
computation section 9 through the set point input 

10 section 1. The controlled variable PV is detected by a 
sensor (not shown) and input to the PID control 
computation section 4 , limit cycle autotuning 
computation section 7, and two-step limit cycle 
computation section 9 through the controlled variable 

15 input section 2. 

When the operator issues an instruction to 
execute limit cycle autotuning, the two-step limit cycle 
computation section 9 is activated to start a series of 
operations (step 101 in Fig. 2). The activated two-step 

20 limit cycle computation section 9 instructs the limit 
cycle autotuning computation section 7 to execute the 
first limit cycle. In response to this instruction, the 
limit cycle autotuning computation section 7 reads out 
the heat -side manipulated variable set point ATh from 

25 the heat-side manipulated variable storage section 5, 
and the cool -side manipulated variable set point ATc 
from the cool- side manipulated variable storage section 
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6 (step 102) . 

The limit cycle autotuning computation section 

7 compares the controlled variable PV with the set point 
SP at the time of execution of the first limit cycle 

5 (step 103). If the controlled variable PV is equal to 
or lower than the set point SP, the limit cycle 
autotuning computation section 7 outputs the heat -side 
manipulated variable set point ATh to the heater (not 
shown) through the manipulated variable output section 3 

10 (step 104). If the controlled variable PV is higher 
than the set point SP, the limit cycle autotuning 
computation section 7 outputs the cool- side manipulated 
variable set point ATc to the cooler (not shown) through 
the manipulated variable output section 3 (step 105). 

15 The limit cycle autotuning computation section 

7 then performs the first ascent/descent extreme value 
detection processing (step 106). Fig. 4 shows the first 
ascent/descent extreme value detection processing by the 
limit cycle autotuning computation section 7. Fig. 5 

20 explains the first ascent /descent extreme value 
detection processing. 

First of all, the limit cycle autotuning 
computation section 7 calculates a deviation Er in the 
current control period on the basis of the set point SP 

25 and controlled variable PV (step 201 in Fig. 4). 

Er = SP - PV ... (25) 

Subsequently, the limit cycle autotuning computation 
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section 7 determines whether or not the following 
inequality holds (step 202): 

I Er I > I Ermax I ... (26) 

where Ermax is the maximum deviation whose initial value 
5 is 0. If inequality (26) holds, the limit cycle 

autotuning computation section 7 sets Ermax = Er, i.e., 
sets the deviation Er in the current control period to 
the maximum deviation Ermax (step 203). 

The limit cycle autotuning computation section 
10 7 then determines according to the following inequality 
whether or not the sign of the deviation Er has changed 
(step 204) . 

ErErO < 0 ... (27) 

where ErO is a. deviation one control period before the 

15 current one. Inequality (27) is used to determine that 
the sign of the deviation Er is reversed, when the 
product of the current deviation Er and the deviation 
ErO one control period before the current one is 
negative. If inequality (27) does not hold, it is 

20 determined that ascent/descent extreme value detection 
is not complete, and the flow returns to step 103. 

When the processing in steps 103 to 105 and 
106 (steps 201 to 204) is repeated for every control 
period, the maximum deviation Ermax is updated with an 

25 increase in the deviation Er. At time tl in Fig. 5, 

inequality (27) holds. When inequality (27) holds, the 
limit cycle autotuning computation section 7 sets Erl = 
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Ermax, i.e., sets the maximum deviation Ermax to a first 
extreme deviation Erl . The limit cycle autotuning 
computation section 7 sets, as a first manipulated 
variable switching elapsed time Thl, the time interval 
5 from the time when inequality (27) has previously held 
to the latest time when the maximum deviation Ermax is 
updated (step 205). Note that when inequality (27) 
holds for the first time, the first manipulated variable 
switching elapsed time Thl is set to 0 . 

10 The limit cycle autotuning computation section 

7 determines whether or not the first ascent/descent 
extreme value detection completion condition has held 
(step 206). In this embodiment, assume that when four 
extreme values of the controlled variable PV are 

15 detected, the first ascent/descent extreme value 

detection completion condition holds. In this case, 
since only one extreme value of the controlled variable 
PV is detected, it is determined that the ascent /descent 
extreme value detection is not complete, the maximum 

20 deviation Ermax is initialized to 0 (step 207), and the 
flow returns to step 103. 

When the processing in steps 103 to 106 is 
repeated for every control period, inequality (27) holds 
again at time t3 in Fig. 5. When inequality (27) holds, 

25 the limit cycle autotuning computation section 7 sets 

Er2 = Erl, Erl = Ermax, and Th2 = Thl, i.e., assigns the 
value of the first extreme deviation Erl to the second 
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extreme deviation Er2, sets the maximum deviation Ermax 
to the new first extreme deviation Erl, and assigns the 
value of the first manipulated variable switching 
elapsed time Thl to the second manipulated variable 
5 switching elapsed time Th2. The limit cycle autotuning 
computation section 7 sets, as the new first manipulated 
variable switching elapsed time Thl, the time interval 
from time tl when inequality (27) has previously held to 
latest time t2 when the maximum deviation Ermax is 

10 updated (step 205). 

The limit cycle autotuning computation section 
7 then determines whether or not the first 
ascent/descent extreme value detection completion 
condition has held (step 206). In this case, since only 

15 two extreme values of the controlled variable PV have 

been detected, it is determined that the ascent /descent 
extreme value detection is not complete, and the maximum 
deviation Ermax is initialized to 0 (step 207). The 
flow returns to step 103. 

20 When the processing in steps 103 to 106 is 

repeated for every control period, inequality (27) holds 
again at time t5 in Fig. 5. When inequality (27) holds, 
the limit cycle autotuning computation section 7 sets 
Er3 = Er2, Er2 = Erl, Erl = Ermax, and Th2 = Thl, i.e., 

25 assigns the value of the second extreme deviation Er2 to 
the third extreme deviation Er3, assigns the value of 
the first extreme deviation Erl to the second extreme 
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deviation Er2, sets the maximum deviation Ermax to the 
new first extreme deviation Erl # and assigns the value 
of the first manipulated variable switching elapsed time 
Thl to the second manipulated variable switching elapsed 
5 time Th2. The limit cycle autotuning computation 
section 7 also sets, as the new first manipulated 
variable switching elapsed time Thl, the time interval 
from time t3 when inequality (27) has previously held to 
latest time t4 when the maximum deviation Ermax is 

10 updated (step 205). 

The limit cycle autotuning computation section 
7 then determines whether or not the first 
ascent/descent extreme value detection completion 
.condition has held (step 206). In this case, since only 

15 three extreme values of the controlled variable PV have 
been detected, it is determined that the ascent/descent 
extreme value detection is not complete, and the maximum 
deviation Ermax is initialized to 0 (step 207). The 
flow returns to step 103. 

20 When the processing in steps 103 to 106 is 

repeated for every control period, inequality (27) holds 
again at time t7 in Fig. 5. When inequality (27) holds, 
the limit cycle autotuning computation section 7 sets 
Er3 = Er2, Er2 = Erl, Erl = Ermax, and Th2 = Thl, and 

25 sets, as the new first manipulated variable switching 
elapsed time Thl, the time interval from time t5 when 
inequality (27) has previously held to latest time t6 
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when the maximum deviation Ermax is updated (step 205). 

The limit cycle autotuning computation section 
7 then determines whether or not the first 
ascent /descent extreme value detection completion 
5 condition has held (step 206). In this case, since four 
extreme values of the controlled variable PV have been 
detected, it is determined that the ascent /descent 
extreme value detection is complete, and the maximum 
deviation Ermax is initialized to 0 (step 208). The 

10 flow advances to step 107. 

As is obvious from Fig. 5, although the number 
of extreme values of the controlled variable PV which 
are required to calculate PID parameters is basically 
three, since there is a possibility that the first 

15 extreme value may be an inappropriate value for 

parameter calculation, four extreme values of the 
controlled variable PV are detected. 

Upon completion of the first ascent/descent 
extreme value detection processing, the limit cycle 

20 autotuning computation section 7 calculates the 

amplitude Al of the controlled variable PV in the first 
limit cycle according to the following equation, and 
stores the resultant data in the first limit cycle 
information storage section 10 (step 107). 

25 Al = I Er2 - Erl I ... (28) 

Fig. 6 shows how the amplitude Al is calculated. The 
first limit cycle is terminated in the above manner. 
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Subsequently, the two-step limit cycle 
computation section 9 instructs the limit cycle 
autotuning computation section 7 to execute the second 
limit cycle. In response to this instruction, the limit 
5 cycle autotuning computation section 7 reads out the 
heat -side manipulated variable set point ATh from the 
heat-side manipulated variable storage section 5, and 
also reads out the cool -side manipulated variable set 
point ATc from the cool -side manipulated variable 

10 storage section 6. In addition, the limit cycle 

autotuning computation section 7 reads out, from the 
manipulated variable change ratio storage section 8 , the 
change instruction information CI for instructing which 
one of the heat -side manipulated variable set point ATh 

15 and the cool-side manipulated variable set point ATc is 
to be changed in the second limit cycle, and the 
manipulated variable change ratio M (step 108 in 
Fig. 3). 

If it is known from information on control 
20 apparatus design that there is a difference between the 
heating capacity on the heat side and the cooling 
capacity on the cool side, the change instruction 
information CI is set in advance before the execution of 
autotuning so as to decrease the manipulated variable 
25 set point on the side with higher capacity in order to 
correct this difference. If, therefore, the heating 
capacity on the heat side is higher, the change 
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instruction information CI instructs to change the 
heat -side manipulated variable set point ATh to //ATh. 
If the cooling capacity on the cool side is higher, the 
change instruction information CI instructs to change 
5 the cool -side manipulated variable set point ATc to 

MATc. The manipulated variable change ratio H is set 
as a coefficient of about 0*6 to 0.8. 

The limit cycle autotuning computation section 
7 then determines on the basis of the change instruction 

10 information CI which one of the cool- side manipulated 
variable set point and the heat -side manipulated 
variable set point is to be changed (step 109). If the 
change instruction information CI instructs to change 
the cool -side manipulated variable set point ATc, the 

15 limit cycle autotuning computation section 7 compares 
the controlled variable PV with the set point SP (step 
110). If the controlled variable PV is equal to or 
lower than the set point SP, the limit cycle autotuning 
computation section 7 outputs the heat -side manipulated 

20 variable set point ATh to the heater (step 111). If the 
controlled variable PV is higher than the set point SP, 
the limit cycle autotuning computation section 7 changes 
the cool -side manipulated variable set point to //ATc 
and outputs it to the cooler (step 112). 

25 The limit cycle autotuning computation section 

7 then performs the second ascent/descent extreme value 
detection processing (step 113). The second 
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ascent /descent extreme value detection processing is 
similar to the first ascent/descent extreme value 
detection processing described with reference to Figs. 4 
and 5. The second ascent/descent extreme value 
5 detection processing differs from the first 

ascent/descent extreme value detection processing only 
in the following point . Of four detected extreme 
values, a deviation in the latest extreme value is set 
as a fourth extreme deviation Er4 (corresponding to Erl 

10 in Fig. 5); a deviation in the second newest extreme 

value, as a fifth extreme deviation Er5 (corresponding 
to Er2 in Fig. 5); a deviation in the third newest 
extreme value, as a sixth extreme deviation Er6 
(corresponding to Er3 in Fig. 5); the time interval from 

15 the time when the sign of the deviation Er is reversed 
immediately before the fourth extreme deviation Er4 to 
the time when the fourth extreme deviation Er4 is 
obtained, as a third manipulated variable switching 
elapsed time Th3 (corresponding to Thl in Fig. 5); and 

20 the time interval from the time when the sign of the 
deviation Er is reversed immediately before the fifth 
extreme deviation Er5 to the time when the fifth extreme 
deviation Er5 is obtained, as a fourth manipulated 
variable switching elapsed time Th4 (corresponding to 

25 Th2 in Fig. 5) . 

Upon completion of the second ascent/descent 
extreme value detection processing, the limit cycle 
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autotuning computation section 7 calculates the 
amplitude A2 of the controlled variable PV in the second 
limit cycle according to the following equation, and 
stores the resultant data in the second limit cycle 
5 information storage section 11 (step 114) 

A2 = I Er5 - Er4 I ... (29) 

The limit cycle autotuning computation section 
7 also stores, in the second limit cycle information 
storage section 11, the manipulated variable switching 

10 elapsed time Th_Heat from the instant at which the 

manipulated variable MV is switched from the cool- side 
manipulated variable set point MATc to the heat -side 
manipulated variable set point ATh in the second limit 
cycle to the instant at which the controlled variable PV 

15 reaches the minimum value, and the manipulated variable 
switching elapsed time Th_Cool from the instant at which 
the manipulated variable MV is switched from the 
heat -side manipulated variable set point ATh to the 
cool- side manipulated variable set point MATc to the 

20 instant at which the controlled variable PV reaches the 
maximum value (step 114). If the fourth extreme 
deviation Er4 is negative (the controlled variable PV is 
the maximum value), the third manipulated variable 
switching elapsed time Th3 is Th_Cool, and the fourth 

25 manipulated variable switching elapsed time Th4 is 
Th_Heat. If the fourth extreme deviation Er4 is 
positive (the controlled variable PV is the minimum 
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value), the third manipulated variable switching elapsed 
time Th3 is Th_Heat, and the fourth manipulated variable 
switching elapsed time Th4 is Th_Cool. 

If it is determined in step 109 that the 
5 change instruction information CI instructs to change 
the heat -side manipulated variable set point ATh, the 
limit cycle autotuning computation section 7 compares 
the controlled variable PV with the set point SP (step 
115). If the controlled variable PV is equal to or 

10 lower than the set point SP, the limit cycle autotuning 
computation section 7 changes the heat -side manipulated 
variable set point to M ATh, and outputs it to the 
heater (step 116). If the controlled variable PV is 
higher than the set point SP, the limit cycle autotuning 

15 computation section 7 outputs the cool -side manipulated 
variable set point ATc to the cooler (step 117). 

Subsequently, the limit cycle autotuning 
computation section 7 performs the second ascent/descent 
extreme value detection processing (step 118), and 

20 stores the amplitude A2 of the controlled variable PV in 
the second limit cycle, the manipulated variable 
switching elapsed time Th_Heat , and the manipulated 
variable switching elapsed time Th_Cool in the second 
limit cycle information storage section 11 (step 119). 

2 5 The processing in steps 118 and 119 is the same as that 
in steps 113 and 114. Fig. 6 shows how the amplitude A2 
is calculated. 
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Upon completion of step 114 or step 119, the 
PID parameter calculation section 12 calculates average 
PID parameters on the basis of the amplitude A2 of the 
controlled variable PV, the manipulated variable 
5 switching elapsed time Th_Heat, and the manipulated 

variable switching elapsed time Th_Cool which are stored 
in the second limit cycle information storage section 11 
(step 120) . 

That is, when the heat -side manipulated 

10 variable set point is changed to MATh in the second 
limit cycle, the PID parameter calculation section 12 
calculates the average proportional band Pb_Ave 
according to equation (15). When the cool-side 
manipulated variable set point is changed to /XATc in 

15 the second limit cycle, the PID parameter calculation 
section 12 calculates the average proportional band 
Pb_Ave according to equation (16). When either the 
heat -side manipulated variable set point is changed to 
U ATh or the cool -side manipulated variable set point is 

20 changed to /XATc in the second limit cycle, the PID 

parameter calculation section 12 calculates the average 
integral time Ti_Ave according to equation (19), and 
calculates the average derivative time Td__Ave according 
to equation (20). 

25 Subsequently, when the heat -side manipulated 

variable set point is changed to tt ATh in the second 
limit cycle, the PID parameter calculation section 12 



- 31 - 



calculates the proportional band ratio k according to 
equation (13). When the cool-side manipulated variable 
set point is changed to /XATc in the second limit cycle, 
the PID parameter calculation section 12 calculates the 
5 proportional band ratio K according to equation (14) 
(step 121) . 

The PID parameter calculation section 12 
calculates PID parameters including the proportional 
band Pb, integral time Ti, and derivative time Td on 

10 each of the heat side and the cool side (step 122). 

When either the heat -side manipulated variable set point 
is changed to /XATh or the cool-side manipulated 
variable set point is changed to MATc in the second 
limit cycle, the PID parameter calculation section 12 

15 calculates the heat-side proportional band Pb_Heat, 
heat -side integral time Ti_Heat, and heat -side 
derivative time Td_Heat according to equations (17), 
(21), and (23), respectively, and calculates the 
cool- side proportional band Pb_Cool, cool -side integral 

20 time Ti_Cool, and cool -side derivative time Td_Cool 
according to equations (18), (22), and (24), 
respectively . 

Finally, the PID parameter calculation section 
12 sets the calculated heat-side proportional band 

25 Pb_Heat, heat-side integral time Ti_Heat, heat-side 
derivative time Td_Heat, cool-side proportional band 
Pb_Cool, cool-side integral time Ti_Cool, and cool-side 
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derivative time Td_Cool in the PID control computation 
section 4 (step 122). With the above operation, the PID 
parameter calculation processing is complete, and limit 
cycle autotuning is terminated. 
5 In normal control operation after limit cycle 

autotuning, the PID control computation section 4 
performs PID control computation to calculate the 
manipulated variable MV for every control period on the 
basis of the set point SP input from the set point input 

10 section 1 and the controlled variable PV input from the 
controlled variable input section 2. At this time, if 
the manipulated variable MV one control period before 
the current one is equal to or higher than a 
predetermined value (e.g., 0%), the PID control 

15 computation section 4 determines the heat mode, and 
calculates the manipulated variable MV to make the 
controlled variable PV agree with the set point SP 
according to the following equation. The PID control 
computation section 4 outputs the calculated manipulated 

20 variable MV to the heater through the manipulated 
variable output section 3 . 

MV = ( £ /Pb_Heat){l + ( l/Ti_Heats ) 

+ Td_Heats} (SP - PV) ...(30) 
where s is a Laplace operator, and the constant C is, 

25 for example, 100. 

If the manipulated variable MV one control 
period before the current one is less than the 
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predetermined value, the PID control computation section 
4 determines the cool mode, and calculates the 
manipulated variable MV according to the following 
equation. The PID control computation section 4 outputs 
5 the calculated manipulated variable MV to the cooler 
through the manipulated variable output section 3. 
MV = ( C /Pb_Cool){l + (l/Ti_Cools) 

+ Td_Cools} (SP - PV) ...(31) 
As described above, according to this 

10 embodiment, at the time of execution of limit cycle 
autotuning, the first limit cycle in which the 
predetermined heat -side manipulated variable set point 
ATh and the predetermined cool- side manipulated variable 
set point ATc are alternately output is generated to 

15 detect the first control response corresponding to the 
first limit cycle, and the second limit cycle in which 
either the heat -side manipulated variable set point ATh 
or the cool-side manipulated variable set point ATc is 
changed on the basis of the predetermined change 

20 instruction information CI and manipulated variable 
change ratio H is generated to detect the second 
control response corresponding to the second limit 
cycle. Control parameters (PID parameters) are then 
calculated on the basis of the first and second control 

25 responses. Even if, therefore, no limit cycle can be 
generated by a heating actuator or cooling actuator 
alone, control parameters suitable for the heat and cool 
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modes can be obtained. 

In this embodiment, the first amplitude of the 
controlled variable is detected as the first control 
response, and the second amplitude of the controlled 
5 variable, a heat-side elapsed time, and a cool-side 
elapsed time are detected are detected as the second 
control response. The ratio between the heat-side 
process gain and the cool- side process gain is obtained 
on the basis of the first and second amplitudes, and a 

10 proportion band in each of the heat mode and the cool 
mode is calculated from this ratio. Integral and 
derivative times common to the heat and cool modes are 
then calculated from the average of the heat -side 
elapsed time and cool -side elapsed time. This makes it 

15 possible to easily obtain control parameters including a 
proportional band, integral time, and derivative time. 
In addition, using the average of a heat -side elapsed 
time and cool -side elapsed time for the calculation of 
an integral time and derivative time can prevent the 

20 occurrence of discontinuous control operation when the 
heat mode is switched to the cool mode or the cool mode 
is switched to the heat mode in normal control operation 
after the completion of autotuning. 
Second Embodiment 

25 The second embodiment of the present invention 

will be described next. The principle of this 
embodiment will be described first. In the first 
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embodiment, the change instruction information CI for 
instructing to change either the heat -side manipulated 
variable set point ATh or the cool -side manipulated 
variable set point ATc in the second limit cycle and the 
5 manipulated variable change ratio fi are set in advance. 
It is difficult for the operator who has no expert 
knowledge about control to make such settings . This 
embodiment therefore provides an arrangement for making 
a heat/cool control apparatus automatically determine 

10 change instruction information CI and a manipulated 
variable change ratio ii . 

In an automatic determination method, a 
heat-side maximum deviation Erh (Erh < 0) and a 
cool -side maximum deviation Ere (Ere > 0) which are 

15 caused in the first limit cycle are detected. The ratio 
between the maximum deviations Erh and Ere reflects the 
balance between the heating capacity on the heat side 
and the cooling capacity on the cool side, and hence the 
change instruction information CI and manipulated 

20 variable change ratio IX may be automatically determined 
on the basis of this ratio. 

More specifically, if I Ere I ^ I Erh I holds, 
the manipulated variable change ratio l± may be 
calculated according to the following equation to 

25 generate the change instruction information CI for 

instructing to change the cool -side manipulated variable 
set point ATc. 
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U = I Erh/Erc I ... (32) 

If I Erh I < I Ere I holds, the manipulated variable 
change ratio U may be calculated according to the 
following equation to generate the change instruction 
5 information CI for instructing to change the heat -side 
manipulated variable set point ATh. 

U = I Erc/Erh I ... (33) 

The arrangement of the heat /cool control 
apparatus according to this embodiment will be described 

10 on the basis of the above principle. Fig. 7 shows the 

arrangement of the heat /cool control apparatus according 
to the second embodiment of the present invention. The 
same reference numerals as in Fig. 1 denote the same 
parts in Fig. 7. A limit cycle autotuning computation 

15 section 7a has the same function as that of the limit 
cycle autotuning computation section 7 in Fig. 1, and 
also has a function of storing the heat -side maximum 
deviation Erh and cool -side maximum deviation Ere in a 
first limit cycle information storage section 10 after 

20 the first limit cycle. In place of the manipulated 

variable change ratio storage section 8 in Fig. 1, this 
embodiment also has an manipulated variable change ratio 
calculation section 8a which determines the change 
instruction information CI and manipulated variable 

25 change ratio U on the basis of the first response 

corresponding to the first limit cycle. The manipulated 
variable change ratio calculation section 8a implements 
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a manipulated variable change ratio calculation means. 

The operation of the heat/cool control 
apparatus in Fig. 7 will be described next with 
reference to Fig. 8. Steps 102 to 105 in Fig. 8 show a 
5 sequence for generating the first limit cycle; steps 106 
and 107a, a sequence for detecting the first control 
response; step 123, a sequence for calculating a 
manipulated variable change ratio; step 108a and steps 
109 to 112, and steps 115 to 117 in Fig. 3, a sequence 

10 for generating the second limit cycle; steps 113, 114, 
118, and 119, a sequence for detecting the second 
control response; and steps 120 to 122, a sequence for 
calculating control parameters. 

When the operator issues an instruction to 

15 execute limit cycle autotuning, a two-step limit cycle 
computation section 9 is activated to start a series of 
operations (step 101 in Fig. 8). The processing in 
steps 102 to 106 is the same as that in the first 
embodiment . 

20 After the completion of the first 

ascent /descent extreme value detection processing, the 
limit cycle autotuning computation section 7a calculates 
an amplitude Al of a controlled variable PV in the first 
limit cycle according to equation (28), and stores the 

25 resultant data in the first limit cycle information 
storage section 10. The limit cycle autotuning 
computation section 7a also stores the heat -side maximum 
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deviation Erh and cool- side maximum deviation Ere in the 
first limit cycle in the first limit cycle information 
storage section 10 (step 107a). If a first extreme 
deviation Erl is negative (the controlled variable PV is 
5 the maximum value), the first extreme deviation Erl is 
Erh, and the second extreme deviation Er2 is Ere, If 
the first extreme deviation Erl is positive (the 
controlled variable PV is the minimum value), the first 
extreme deviation Erl is Ere, and the second extreme 

10 deviation Er2 is Erh. 

The manipulated variable change ratio 
calculation section 8a compares the heat -side maximum 
deviation Erh and cool- side maximum deviation Ere stored 
in the first limit cycle information storage section 10. 

15 If I Ere I ^ I Erh I holds, the manipulated variable 
change ratio calculation section 8a calculates the 
manipulated variable change ratio U according to 
equation (32) to generate the change instruction 
information CI for instructing to change the cool-side 

20 manipulated variable set point ATc. If I Ere I < I Erh I 
holds, the manipulated variable change ratio calculation 
section 8a calculates the manipulated variable change 
ratio U according to equation (33) to generate the 
change instruction information CI for instructing to 

25 change the heat -side manipulated variable set point ATh 
(step 123). The manipulated variable change ratio 
calculation section 8a then stores the calculated 
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manipulated variable change ratio U and the generated 
change instruction information CI . 

After the completion of the processing in step 
123, the two-step limit cycle computation section 9 
5 instructs the limit cycle autotuning computation section 
7a to execute the second limit cycle. In response to 
this instruction, the limit cycle autotuning computation 
section 7a reads out the heat -side manipulated variable 
set point ATh from a heat -side manipulated variable 

10 storage section 5 and the cool- side manipulated variable 
set point ATc from a cool-side manipulated variable 
storage section 6. The limit cycle autotuning 
computation section 7a further reads out the change 
instruction information CI and manipulated variable 

15 change ratio ix from the manipulated variable change 
ratio calculation section 8a (step 108a). 

Since the processing in steps 109 to 122 is 
the same as that in the first embodiment, an 
illustration and description thereof will be omitted. 

20 As described above, according to this 

embodiment, the change instruction information CI and 
manipulated variable change ratio U can be 
automatically determined on the basis of the first 
control response. As a consequence, according to this 

25 embodiment, the same effect as that of the first 

embodiment can be obtained. In addition, since there is 
no need to set the change instruction information CI and 
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manipulated variable change ratio U in advance, a 
heat /cool control apparatus can be realized, which 
allows even an operator who has no expert knowledge 
about control to easily use. 
5 In addition, in this embodiment, by detecting 

a heat -side maximum deviation and cool- side maximum 
deviation as the first control response, change 
instruction information and a manipulated variable 
change ratio can be easily determined on the basis of 

10 the heat -side maximum deviation and cool -side maximum 
deviation. In addition, in the control parameter 
calculation sequence, calculating heat-mode integral and 
derivative times from a heat-side elapsed time, and 
cool -mode integral and derivative times from a cool- side 

15 elapsed time makes it possible to obtain more preferable 
integral and derivative times for each of the heat mode 
and the cool mode, thereby obtaining better control 
characteristics . 
Third Embodiment 

20 The third embodiment of the present invention 

will be described next. The principle of this 
embodiment will be described first. In the second 
embodiment , the manipulated variable change ratio tt is 
automatically determined. According to this automatic 

25 determination method, since the difference between the 
heat -side maximum deviation Erh and the cool- side 
maximum deviation Ere in the second limit cycle becomes 



- 41 - 



smaller than the first limit cycle, more ideal 
autotuning can be realized. 

In the first and second embodiments, the 
average integral time Ti_Ave and average derivative time 
5 Td_Ave are used as the integral time Ti and derivative 
time Td common to the heat side and cool side on the 
basis of the premise that the integral time Ti and 
derivative time Td as PID parameters remain unchanged 
both in the heat mode and the cool mode. The adjustment 

10 method based on the average integral time Ti_Ave and 

average derivative time Td_Ave has a merit that when the 
heat mode is switched to the cool mode or the cool mode 
is switched to the heat mode in normal control operation 
after autotuning, the manipulated variable change is 

15 small and no discontinuous control operation occurs. 

Depending on the characteristics of a control 
target, however, the integral time Ti and derivative 
time Td may be preferably set to different values in the 
heat mode and cool mode. This embodiment provides a 

20 heat/cool control apparatus which can obtain better 

control characteristics by obtaining the integral time 
Ti and derivative time Td suitable for the heat mode and 
cool mode, respectively. 

If a manipulated variable change ratio U is 

25 automatically determined as in the second embodiment, a 
response waveform ideal for autotuning appears in the 
second limit cycle. In this case, a manipulated 
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variable switching elapsed time Th_Cool from the instant 
at which a manipulated variable MV is switched from a 
heat -side manipulated variable set point /ZATh (or ATh) 
to a cool- side manipulated variable set point ATc (or 
5 MATc) to the instant at which a controlled variable PV 
reaches the maximum value is almost equal to a cool -side 
process dead time Lp_Cool set when the cool- side 
manipulated variable set point ATc ( //ATc) is output. 
In addition, a manipulated variable switching elapsed 

10 time Th_Heat from the instant at which the manipulated 
variable MV is switched from a cool- side manipulated 
variable set point ATc ( M ATh) to a heat-side 
manipulated variable set point /x ATh (or ATh) to the 
instant at which the controlled variable PV reaches the 

15 minimum value is almost equal to a heat -side process 
dead time Lp_Heat set when the heat -side manipulated 
variable set point ATh ( //ATh) is output. 



switching elapsed times Th__Heat and Th_Cool and 
20 equations (8) and (9), therefore, a heat -side integral 

time Ti__Heat, cool-side integral time Ti_Cool, heat-side 
derivative time Td_Heat, and cool-side derivative time 
Td_Cool can be calculated as follows: 



On the basis of the manipulated variable 



Ti_Heat = 3 Lp_Heat = 0 Th_Heat 



. . .(34) 



25 



Ti_Cool = /3Lp_Cool = j3Th_Cool 



. . . (35) 



Td_Heat = 7Lp_Heat = 7Th_Heat 



. . . (36) 



Td_Cool = 7Lp_Cool = 7Th_Cool 



. . .(37) 
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The arrangement and operation of the heat /cool 
control apparatus of this embodiment are almost the same 
as those of the second embodiment, and differ only in 
the following points. First, there is no need to 
5 calculate the average integral time Ti_Ave and average 
derivative time Td_Ave in step 120 in Fig. 3. Second, 
when calculating PID parameters in step 122, a PID 
parameter calculation section 12 calculates the 
heat -side integral time Ti_Heat and heat -side derivative 

10 time Td_Heat according to equations (34) and (36), 

respectively, and calculates the cool- side integral time 
Ti_Cool and cool- side derivative time Td_Cool according 
to equations (35) and (37), respectively. 

As has been described above, according to this 

15 embodiment, more preferable integral time Ti and 

derivative time Td can be obtained in each of the heat 
mode and the cool mode. As a consequence, in this 
embodiment, the same effect as that in the second 
embodiment can be obtained, and control characteristics 

20 better than those in the first and second embodiments 
can be obtained. 
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